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Abstract 

Both the trapping geometry and the interatomic interaction strength of a dilute ultracold 
fermionic gas can be well controlled experimentally. When the interactions are tuned to strong 
attraction, Cooper pairing of neutral atoms takes place and a BCS superfluid is created. Alter- 
natively, the presence of Feshbach resonances in the interatomic scattering allows populating a 
molecular (bound) state. These molecules are more tightly bound than the Cooper pairs and can 
form a Bose-Einstein condensate (BEC). In this contribution, we describe both the BCS and BEC 
regimes, and the crossover, from a functional integral point of view. In this description, the prop- 
erties of the superfluid (such as vortices and Josephson tunneling) can be derived and followed as 
the system is tuned from BCS the BEC. In particular, we present results for the critical current of 
the superfluid through an optical lattice and link these results to recent experiments with atomic 
bosons in optical lattices. 



PACS numbers: 03.75.-b, 03.75.Lm, 03.75. Ss 



I. THE ULTRACOLD DILUTE FERMI GAS 



When a dilute Bose gas is cooled below the degeneracy temperature, the bosonic atoms all 
condense in the same one-particle state and a Bose-Einstein condensate forms. This has been 
convincingly demonstrated with magnetically trapped, evaporatively cooled atomic gases for 
a multitude of atom species. Moreover, magnetic or optical traps can be equally well loaded 
with fermionic isotopes, such as 6 Li or 40 K. These atoms do not undergo Bose-Einstein 
condensation, but fill up a Fermi sea, as has been demonstrated through the observation of 
the Pauli blocking effect £0] and through a measurement of the total energy of the Fermi 
gas JJ. Very soon after the observation of a degenerate Fermi sea of atoms, researchers 
embarked upon the quest to achieve Cooper pairing in the dilute Fermi gas. Indeed, for 
metals we know that the Fermi sea is unstable with respect to Cooper pair formation. So, 
if the (neutral) atoms in the dilute gas attract each other, a similar instability towards a 
paired state is to be expected. 

The interatomic interactions in ultracold gases are remarkable for two reasons. Firstly, 
the collisions between the atoms can be satisfactorily characterized by a single number, the 
s-wave scattering length a s . For low-energy collisions, the effective interaction potential 
between atoms becomes a contact potential, V{r — r') = g5(r — r'), where g = ATrha s /m with 
m the mass of the atoms. The scattering length can be both positive (leading to interatomic 
repulsion) or negative (attraction). 

Secondly, this scattering length can be tuned by an external magnetic field when a Fesh- 
bach resonance is present |3j . This resonance occurs when the energy of a bound (molecular) 
state in a closed scattering channel becomes equal to the energy of the colliding atoms in 
the open scattering channel. The different channels correspond here to different hyperfine 
states of the trapped atoms, and the distance in energy between these states can be tuned 
with a magnetic field. 

In what follows, we will consider a trapped mixture of 40 K atoms in the 1 9/2, —7/2) and 
1 9/2, — 9/2) hyperfine states. This potassium isotope is fermionic, and the trapped states 
display a Feshbach resonance at B = 202.1 Gauss. When the scattering length is tuned to a 
negative value, the atoms attract and Cooper pairs can form leading to a BCS regime. The 
critical temperature for Cooper pairing can be raised by making the scattering length more 
strongly negative. When the scattering length is large and positive, the molecular state in 
the closed channel is populated, and molecules appear that can be Bose-Einstein condensed 
(the BEC regime). The adaptability of the scattering length allows bringing the gas from 
the BCS regime into the BEC regime or vice versa, and allows studying the interesting 
intermediate 'crossover' regime. 

The first experimental realization of superfluidity of a Fermi gas in the molecular BEC 
regime came in 2003 . A condensate of molecules was convincingly observed. The detec- 
tion of superfluidity in the BCS regime however is much more subtle. In an initial experiment 
0, the superfluid behavior was derived from the hydrodynamic nature of the expansion 
of the cloud, as compared to a ballistic expansion expected for a non-superfluid weakly- 
interacting Fermi gas 0. However, this did not constitute unambiguous proof, since the 
Fermi gas was in the strongly interacting regime. Subsequent experiments probed superflu- 
idity by mapping the pair density onto a molecular condensate density 0] or by spectroscop- 
ically measuring the gap |8| . Yet although these experimental methods clearly demonstrate 
pairing, they do not unambiguously demonstrate superfluid behavior. 
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FIG. 1: Two counterpropagating laser beams form a periodic potential for the atoms. Such optical 
lattices can be loaded with quantum gases, forming a stack of quasi-two dimensional clouds. 



The very recent observation of a lattice of quantized vortices in resonant Fermi gases 
constitutes the first clear demonstration of superfluidity in the BEC/BCS regime. Obser- 
vation of these vortices well in the BCS regime may be difficult since the fermionic density 
penetrates in the core of the vortex in the BCS regime, leading to a loss of contrast in direct 
imaging (HUH U|- Another possibility to demonstrate superconductivity is though the ob 



servation of the Josephson effect l]l in optical lattices. These optical lattices are periodic 



potentials formed by two counterpropagating laser beams, for example in the ^-direction: 

V opt {z) = sE R sm 2 (2irz/\), (1) 

where A is the laser wave length, Er = h 2 / (2mA 2 ) is the recoil energy, and s is the laser 
intensity expressed in units of the recoil energy. Typically, s = 1 — 20, A = 795 nm. The 
atoms collect in the valleys of the optical lattice and form a "stack of pancakes", illustrated 
in Fig. ^ Typically, there are on the order of a few 100 'pancakes' with on the order of 1000 
atoms each. When a superfluid is loaded in such an optical lattice, the system corresponds to 
an array of Josephson junctions. In such an array, the superfluid gas can propagate whereas 
the normal state gas is pinned. This has already been demonstrated for bosonic atoms , 
and has been predicted theoretically for fermionic atoms (l2j, . 

In this contribution, we derive and discuss the critical Josephson current for the flow of the 
superfluid component through an optical lattice. For this purpose, we base ourselves on the 
path-integral description as applied by Randeria and co- workers [H Q to the BEC/BCS 



crossover model of high-T c superconductors. In section II, we give an overview of the 
application of path-integrals to the system of ultracold fermions, and in section III we 
present our results for the critical current. 



II. PATH-INTEGRAL TREATMENT OF THE BEC/BCS CROSSOVER 

The partition function for the atomic Fermi gas is given by the functional integral 

Z = J PVw£>Vw exp{-S/h} (2) 



with an action 



ft/3 





H/3 



+ f dT f ^ fl , ^X,T,TV ; x,T,|V ; x,T,|V ; x,T,T- (3) 



The fermionic fields ip X}T , t/> xt are Grassman variables. The interaction potential, as discussed 
in the previous section, is a contact potential with experimentally adjustable strength g. 
The two hyperfine states are denoted by o =f , j. The functional integral over the Grassman 
variables can be performed analytically only for an action that is quadratic in ■0 X( t)V'x i t- 
In order to get rid of the quartic term in (JHJ) we perform a Hubbard-Stratonovic (HS) 
transformation, introducing auxiliary bosonic fields A x >T and A x T : 

Zcx^ Thk^Ttyw J PA XiT DA x , T exp{-<S//i} (4) 



with 



ft/3 





ft/3 



- Jdr J dx ^A XiT ^ x>Tii ^ XjTiT + A XjT '0 x ,r ) TV'x,T,| + Ax ' r ^ Ax ' T ^ . (5) 



Indeed, performing the functional integral over the HS fields A XjT ,A X)T in (JHJ) brings us back 
to Our goal is an investigation of the superfluid properties of the ultracold Fermi system. 
For a straightforward hydrodynamic interpretation of the Hubbard-Stratonovic fields, it is 
advantageous to work with |A xr | and # XjT . These are related to the original HS field by 
A XiT = | A XiT | exp(i# X;T ). We have restricted the functional integral to A X T = (A X)T )* without 
neglecting any field configurations of importance to the final result. The hydrodynamic 
interpretation of |A X)T | 2 is the density of fermion pairs, whereas ftV x # Xir /m = v XiT can 
be interpreted as the superfluid velocity field. Performing this change of variables in the 
functional integral yields 

Z oc j PVv )(J £>Vv,* J V |A x ,r| V6 X>T exp {-S/h} , (6) 



with 



ft/3 

f f - ( d h 2 1 
S = dr dx ^x,r, CT ( h— - — V x - -v XjT ■ «W X - \i 

o 

ihdOxr h. m x 1 o \ , 

' j ^ V ' AX ' T ' ^ X ' T '^ X ' T ' T + ' AX ' T ' ^ X <^X,T,| + j • ( 7 ) 
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De Palo et al. 17] suggest at this point to introduce additional collective quantum vari- 
ables to extract the fermionic density. However, care must be taken, since when additional 
collective quantum fields are present the problem of double-counting poses itself [l8j], and 
variational perturbation theory has to be applied to avoid double-counting 1^1 • However, in 



the present case it is not necessary to explicitly introduce the additional collective variables 
to obtain information about the atomic density profile In © the integration over the 
fermionic variables can be taken, leading to 

Z cc J V |A XjT | V9^ T exp {-S cS /h} (8) 



with an effective action 

H/3 



S cS = -hti 



In 



-G- 



h 



IA 



1 2 







dr I d-K (9) 



where the inverse propagator can be written as the sum of an inverse 'free fermion propa- 
gator' and a term arising from the superfluidity: 



i — l _ /pi — l 



(10) 



The inverse free fermion propagator is 



^ = Mft!:Ws(— (11) 



and the superfluid part of the propagator can be written as 

§ = o- (~^ Vx >^ ' ^V X J - ax{h |A X)T |) 

/ ihd6y- T 1 „ 1 9 \ , . 

+ ^3 ( y -g- ~ ^V x • v x , r + -m <T \ . (12) 

In these expressions, ao...a^ are the Pauli matrices. Note that if we have an external potential 
V ext {x) present, for example the optical potential or the magnetic trap, this appears in — Gq 1 
as an extra term +a s V ex t(x). The effective action (jHJ) depends on the fields |A XjT | , X)T and 
p X T , C x T . For the former, a saddle point approximation is usually made. For example, a 
good saddle point form when no vortex is present is [l5L : 



A X)T | = A 



constant. 



(SP1) 



The value of the constant for the phase is irrelevant, and the value of A can be extracted by 
extremizing the effective action 8S e s/SA = 0. This yields the well-known gap equation in the 
case of neutral atoms interacting through a contact potential. Alternatively, we proposed in 
Ref. jll| to use a different saddle point approximation to investigate the case of a fermionic 
superfluid containing a vortex parallel to the z-axis: 
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Here, is the angle around the z-axis, and r is the distance to the z-axis. Again, a gap 
equation can be derived for A r by extremizing the action - this gap equation yields a gap 
that depends on the distance to the vortex line (the z-axis). Fixing the total number of 
fermions yields a number equation in which the local density of fermions can be identified 
st r aight forwar dly 

Consider first the simplest saddle point approximation, (jSPlj) . The saddle point result 
for the action in this case is 



|A| 2 f dk 



SP1 " 9 J (2tt)» 



In 




(13) 



Two unknowns are the chemical potential fi and the value of constant A, the gap. The 
chemical potential is obtained by fixing the particle density In the BCS limit, fi — > Ep 
whereas in the BEC limit, the chemical potential goes to the binding energy of the molecule, 
— > h 2 /(ma 2 ). In the intermediate regime, there is a smooth crossover between the two 
limiting values. The gap A is found by extremizing the saddle point action, 5S sp i/5A = 0. 
The result is shown for different temperatures in figure [21 At temperature zero, the gap 
depends exponentially on the scattering length as we expect from the BCS theory. As 
the temperature is raised, the gap decreases, reaching zero at a certain temperature. In 
the BCS limit, the superfluidity is destroyed by breaking up Cooper pairs, so the critical 
temperature corresponds to the temperature where A = 0. However, in the BEC limit, 
superfluidity is destroyed through phase fluctuations, and one cannot extract the critical 
temperature from the results shown in figure El It becomes necessary to include fluctuations 
around the saddle point value (jSPl|) and expand the effective action up to second order 
in these fluctuations around the saddle point value. This second-order expansion yields an 
action that is quadratic in the fluctuation variables and that can be integrated analytically. 
For fluctuations around the saddle point (|SP1|) this was done by Randeria and co-workers, 
who obtained a corrected value of the critical temperature that in the BEC limit becomes 
independent of l/(A;^a s ). More recently, the effects of fluctuations in the superfluid regime, 
in the context of a diagrammatic expansion of the thermodynamic potential in refs. (21. 22l|. 



III. DETERMINATION OF THE CRITICAL CURRENT 
A. Effective action in the optical potential 

The path-integral method outlined in the previous section has been applied before to 
describe vortices in a superfluid Fermi gas ll|| and to describe the propagation of a superfluid 



Fermi gas in an optical potential 12]. When an optical lattice (JTJ) is present along the in- 
direction, we can decouple the (free) motion in the x, y-plane from the (tunneling) motion 
in the z-direction. To make this decoupling clear in the notations, we write the partition 
function of the system as 

Z = J ViP jja {r)ViP j>a {r) exp {-S/h} , (14) 
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FIG. 2: The saddle point value A is shown as a function of the interaction strength parameter 
l/(fci?o«) f° r different values of the temperature. In the BCS regime this corresponds to the BCS 
gap and it vanishes at the critical temperature. The inset shows the temperature at which A = 0. 
In the BEC regime, fluctuations around the saddle point need to be taken into account to obtain 
the correct critical temperature 
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with r = {x, y, r} and 



5 = E^ + E 5 . 



tunnel 



(15) 



The action functional for the gas in layer j separately is given by 

h/3 

Sj = j dT j J2a^jA r ) [fr-^ ~ ^~ V x _ V ext(j) - fl \ ll>j,a( 



hft 



+ / dT \ dx ^i,T(^)^,i( r )^,i( r )^,T( r )- 



(16) 



This is the two-dimensional version of the action functional (0), supplemented with a layer 
index. Moreover, there is an external potential V ext (j) acting on each layer. This can a 
parabolic potential in addition to the optical potential itself. The tunneling of atoms from 
one layer to another is described by 



hf3 

s tujanei x = J dr J d ^ ff (r)^ +lj(T (r) + ^ i+V (r)^(r)] , (17) 





where the tunneling ener gy U to bring an atom from one well of the optical potential to the 
next was derived in Ref. 12311 : 



h = sE 



In 2 



R 



- 1 



-V5(tt/2) 2 



(18) 
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For this particular decomposition of the action functional in intralayer contributions and 
tunneling contributions, we can perform the same analysis as described in the previous 
section. A Hubbard-Stratonovich transformation gets rid of the four-operator term and 
introduces the HS fields \Aj\,6j, after which the integration over fermionic variables is 
performed. The final result for the effective action can again be written as the sum of 
contributions independent of t 1 and tunneling contributions: 



S, 



eff 



The effective action for layer j is 



effj 



Ectunnel 

j 



(19) 



>cffj 



— htr 



In 



i\ i 



sp 



h 



dr I dx 



I A, I 



(20) 



with 



O"0 



h 



Wt 



°3 



2m 



Vl-V ext {j)-ii )-a l (h\A j \). 



(21) 



The tunneling contributions in the effective action can be treated perturbatively. In that 
framework, the saddle-point values |Aj| can be extracted from the gap equation of each 
layer separately, and the chemical potential p is obtained from the number equation. In 
each layer j, there is an 'effective' chemical potential V ext (j) + p fixing the local density pj 
in layer j. Based on these results for the layers, the lowest-order perturbative expansion of 
the action with respect to the tunneling part (t x ) yields 



otunnel 
°eff, j-^j+l 



hp 

-/*■ 





with, 



T, 



dx T. 



3-+3+1 



cos[^ 



i+i 



(22) 



(23) 



2nh 2 pj/m + Eb 

where E^ is the binding energy of the molecule. This molecular binding energy appears 
through the gap equations and can be derived from scattering theory in reduced dimension- 
ality. It is given by [24^ : 



Ef> = 0.583^E R exp 



A 



2tt s- 



-1/4, 



(24) 



It is important to note that the binding energy depends on the intensity and wavelength of 
the lasers generating the optical potential. More intense laser beams or smaller wavelengths 
confine the gas more strongly in the optical lattice and alter the binding energy of the 
resonant molecules. A more detailed determination of the molecular binding energy in an 
optical lattice, taking into account molecules formed from atoms in neighboring lattice sites, 
is given in Ref. IIF 
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B. Coupled density-phase equations 



The equations of motion for the remaining variables (density pj and phase 9j in layer 
j) can be derived from the effective action (|19 )) -([22 )) through the extremum conditions 
5S e ff/59j = and the number equation. This leads to the equations reported by the present 
authors and M. Wouters in Ref. Ilj 



4m 

+ Tjj_x sin (9j - 9j-i) 

~ Tj+ij sin (9 j+ i - 9j) , (25) 



and 



i2 



-4 = ? L + ^(J')-/' 

2 8m 

--^±008(9,-9^). (26) 

In these equations, we have introduced the possibility of applying an external potential 
Vext(j) varying over the layers. Here, we investigate the case with a constant phase difference 
9j+i — 9j = A9 and a smoothly varying density pj + \ « pj. This situation corresponds to a 
uniform flow of superfluid through the lattice. Then equations (|25|l . (p6*jl simplify to 

^ = (27) 



In the BEC case, Ef, ^> h 2 pj/m and we retrieve the equations describing a conventional 
Josephson junction array. However, on the BCS side, the tunneling coefficients start to 
depend on pj, as Eb and h 2 pj/m become comparable. 



C. Critical Josephson current and critical velocity 

Equation ()27|) states that the current density J is proportional to sin(A6 ) ). This is similar 
to the first Josephson equation, 

J=J c sin(A0). 

The second equation, ()28j) . fixes the time derivative of the phase difference to the differences 
in potential energy of the gas in the different wells. Near A# = 0,7r/2 equation (|28|) yields 
the second Josephson equation —hd(A9)/dt = AU where AU is the difference in energy 
between consecutive wells for the molecules. From the first Josephson equation we can 
extract a critical current by setting A# = n/2, as has been done for bosonic atoms in the 
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FIG. 3: The critical Josephson velocity through the optical lattice, in microns/ms, is shown as a 
function of the interaction strength parameter l/(kpa s ) and the Fermi wavevector of the 2D fermi 
gas in the valleys. In the upper half (a s > 0) of the figure the BEC regime develops, and in the 
lower half (a s < 0) the BCS regime appears. The optical lattice parameters are A = 795 nm and 
s = 3. The atom parameters correspond to 40 K atoms in the 1 9/2, —7/2) and 1 9/2, —9/2) hyperfine 
states, having a Feshbach resonance at B = 202.1 Gauss. 



optical lattice in Ref. [26]. This yields a critical current density for Josephson tunneling 
from layer to layer: 

2 t 2 

Jc = h2nh 2 p,/m + E b Pj ' (29) 

The layers are separated by a distance A/2. From J c we can then extract the critical velocity 
for the fermionic atoms through the optical lattice, 

tl " (30) 



2iTh 2 pj/m + E b H 

The critical velocity of the fermionic superfluid depends on the scattering length a s , via 
the binding energy of the Feshbach resonant molecule, E b . The critical velocity also depends 
on the density (or, equivalently, the Fermi wave vector). In Figure El we show the results for 
the critical velocity (expressed in microns per millisecond), as a function of kp and of the 
interaction parameter l/(kpa s ). In the region l/(kpa s ) > we are in the molecular BEC 
regime, and E b ^> h 2 pj/m. The critical velocity in the BEC regime is roughly proportional 
to t\/E b . In the region l/{kpa s ) < 0, the BCS regime of Cooper pairs arises, and the 
result for the critical velocity becomes nontrivial. For each fixed value of l/(kpa s ) < 0, 
there appears a maximum as a function of kp. This maximum occurs when E b « h 2 pj/m, 
minimizing the denominator in (J3Uj) . 

Although superfluid gases of bosonic atoms have already been studied in optical lattices 
13i 2?J , superfluid Fermi gases have to this moment not been loaded in optical lattices. Also 
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no molecular condensates have been placed in optical lattices. For atomic condensates, a 
critical velocity could be determined experimentally |22j], and was found to vary between 
0.2 and 1.2 /iin/ms for 87 Rb atoms. This is comparable to the velocities that we predict for 
(fermionic) 40 K in the same A = 795 optical potential. Thus, the superffuid regime of paired 
fermionic atoms in an optical lattice should be accessible experimentally. 

IV. CONCLUSIONS 

The path-integral description of ultracold fermionic atoms interacting through a tunable 
contact potential allows to describe vortex configurations and other non-ground state config- 
urations through a judicious choice of saddle point. We apply this formalism to the case of 
a fermionic gas in an optical potential. When the fermionic gas is in the superfluid regime, 
the layers of gas in the optical potential form a Josephson junction array. Equations of 
motion for the density and phase in each layer are obtained and applied to the case where 
the phase difference between consecutive layers is constant. This permits the derivation of a 
critical velocity of the superfluid flow through the optical potential. Although these results 
are strictly speaking derived for T = 0, in the experiments the temperature can typically be 
brought down well below the degeneracy temperature so that we believe our results will be 
relevant to the experiments with optical lattices. 
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